Global Control Methods for GHZ State Generation on 1-D Ising Chain 
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We discuss how to prepare an Ising chain in a GHZ state using a single global control field only. 
This model does not require the spins to be individually addressable and is applicable to quantum 
systems such as cold atoms in optical lattices, some liquid- or solid-state NMR experiments, and 
many nano-scale quantum structures. We show that GHZ states can always be reached asymptoti- 
cally from certain easy-to-prepare initial states using adiabatic passage, and under certain conditions 
finite-time reachability can be ensured. To provide a reference useful for future experimental imple- 
mentations three different control strategies to achieve the objective, adiabatic passage, Lyapunov 
control and optimal control are compared, and their advantages and disadvantages discussed, in par- 
ticular in the presence of realistic imperfections such as imperfect initial state preparation, system 
inhomogeneity and dephasing. 
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I. INTRODUCTION 

Spin chains are an important theoretical model to un- 
derstand properties of many-body systems such as quan- 
tum phase transitions pQ, and have become popular as 
a possible model for quantum computation (QC) and 
communication 2 . Various types of spin-spin inter- 
actions including the Heisenberg, XY and Ising model 
have been discussed both analytically and numerically 
Among these, the Ising model is one of the most ubiq- 
uitous, arising in many different settings from atoms in 
optical lattices [21 E] , to NMR systems [5] to ion traps [B] 
and polar molecules [7]. For many of these systems ad- 
dressing individual spins selectively is extremely difficult, 
limiting the type of control we can implement. For in- 
stance, for cold atoms in an optical lattice, addressing 
individual atoms with an external laser is very difficult 
as the waist of the controlling laser is on the scale of 
many lattice sites. One way to circumvent this problem 
is by eliminating the need for local addressing, i.e., by 
using only control fields that act globally on all the spins 
at once. 

One type of global control was first proposed in [SI [5] , 
where it was demonstrated that universal QC can be 
achieved by controlling the qubits collectively, and since 
then extensive work has been done on this approach |10l - 
IT5] . However, in all of these proposals, the spin- 
spin interaction Hamiltonian and the globally controlled 
Hamiltonian are not sufficient to realize universal QC, 
and additional resources are necessary. For example, 
in the spin chain consists of two types of qubits, A 
and B, arranged in an alternating, repeating pattern: 
ABABAB • • • , and it is assumed that the two types of 
spin-spin interactions Hab and Hba can be switched on 
and off as needed, which is a demanding experimental 
requirement. In an improved scheme [10j . the coupling 
Hamiltonian can be kept unchanged but the transition 



frequencies of each qubit must be tuned individually, 
which is still too difficult to implement for cold atoms 
in optical lattices at this time. In another scheme [11) . 
the need for such "individual tuning ability" is avoided 
but at the expense of requiring a chain with a repeat- 
ing pattern of four types of qubits, raising a high de- 
mand for the preparation process. If the goal is universal 
quantum computation then such extra requirements are 
necessary as a controllability argument shows that the 
system Hamiltonian together with only a global control 
Hamiltonian do not generate the full Lie algebra su(2 JV ) 
but only a proper subalgebra, rendering the system un- 
controllable. One interesting question therefore is what 
further assumptions are necessary for universal QC, a 
question that has been addressed in various recent pa- 
pers [51411] . However, there are many tasks that do not 
require controllability, and we can ask what interesting 
tasks we can perform using global control only, without 
additional resources. This is the focus of this article. 

In particular we demonstrate that global control alone 
is sufficient to steer an Ising chain from a certain ini- 
tial product state to a GHZ state [2]. Such GHZ states 
(or multi-qubit Cat states or NOON states, as they are 
variously called) are of utmost importance for improv- 
ing frequency standards beyond the classical realm [T5] . 
and could result in highly sensitive magnetometers |16| . 
This is why recently there has been extensive interest in 
their experimental realization, for example, in ion traps 
|17) and with multiple nuclear spins in a molecule |18j . 
However, a fast/ non-adiabatic (so as to be robust to 
decoherence) method of generating these systems with 
minimal control, such as global fields on an Ising chain, is 
still an open problem, and will be a very important mile- 
stone for realizing quantum enhanced sensing and stan- 
dards. We again emphasize that all control schemes rely- 
ing purely on global control are useful for experiments on 
the systems where individual addressability is not avail- 
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able. We show that with global control, certain GHZ 
states can be reached in finite time from a given, easy- 
to-prepare product state, i.e., that there always exists 
a control that achieves the task, and consider and com- 
pare three different methods to design suitable controls: 
adiabatic passage, Lyapunov control design and optimal 
control. One of the most interesting issue of our proposal 
is that in all three methods the initial and final states are 
the eigenstates of the Hamiltonian and hence the output 
state does not evolve anymore. This makes it convenient 
since no fine-tuning of the control time is necessary and 
the output state can be saved for further tasks. 

The article is organized as follows: In Sec. II the model 
and control problems are defined. In Sec. Ill, controllabil- 
ity of the system, or rather lack of it, and reachability are 
discussed. In Sec. IV three methods for GHZ state gener- 
ation are considered in detail, namely, adiabatic passage, 
Lyapunov control and optimal control. Finally, the effect 
of various procedural imperfections such as imperfect ini- 
tialization, inhomogeneity and decoherence are studied in 
Sec. V, followed by a brief summary and discussion of 
the results in Sec VI. 



II. ISING MODEL AND CONTROL PROBLEM 

In the following we consider a 1-D spin chain of length 
N with the Hamiltonian 

H f (t) = J[H + f(t)H 1 ], (1) 

where, Hq models the fixed interaction between neigh- 
boring spins with a fixed strength J, and Hi the control 
interaction, corresponding to an applied external field 
f(t). We choose Hq and Hi such that the total Hamilto- 
nian is a uniform nearest-neighbor Ising interaction 
in a transverse time-dependent magnetic field 

N-l N 

Hq = ^ Z n Z n+ i, Hi = X n , (2) 

n— 1 n— 1 



GHZ states 

l^ 1} ) = ^(|0---0) + 11-.-1)) (4a) 

l^ 2) ) = ^(|0---0)-|l---l)) (4b) 

lV>f> = 75(10101 •■■) + |1010 •••)) (4c) 

|V>f> = ^(10101 •■■)- |1010 ■■•)), (4d) 

starting in the ground state \tpf/°^} °f the Hamiltonian 
Hf(0), i.e., our aim is to find a magnetic field f(t) such 
that the system states evolves into one of the entan- 
gled GHZ states given in Eq. Q under the action of 
the resulting Hamiltonian, starting from the initial state 

i4 /o) >- 

In the absence of the magnetic field, /o = 0, the 
ground state of the Hf(0) is two- fold degenerate, spanned 

by {l^lVf}} for J < 0, and {|^ 3) ), |^ 4) )} for 
J > 0. Although the target states are ground states 
of the Hamiltonian Hq, due to the degeneracy cooling 
alone does not suffice to prepare the system in any of 
the states Q. Rather, simple cooling will result in the 
system being left in a mixture of different ground states, 
which is not useful. On the other hand, in the presence 
of the magnetic field the ground state of the total Hamil- 
tonian Hf(0) with /o ^ is non-degenerate, and if the 
system is cooled in the presence of a global field fo along 
the x-axis, it will be initialized in the ground state |V ; c/°' > ) 
of Hf(0). As |/o| — > oo we have 

lim |V^ /o) ) = | + ...+), (5a) 

J/o->--oo 

km \4 M ) = | -...-), (5b) 

J/o-H-oo 

where |±) = ^(|0) ± |1)), and |0) and |1) are the eigen- 
states of Z, i.e., for sufficiently large |/o| we can assume 

IV'o ) to be approximately equal to \i/Jq) = | + . . . +) 
or \ipQ~) = | — ...—). Choosing |/ | large also ensures a 
sufficiently large energy gap for efficient cooling. 



where V, V, Z are the Pauli matrices. Practically, one 
can think of f(t) as a time-dependent global magnetic 
field in the cc-direction that causes all spins to rotate si- 
multaneously, while all spins are constantly coupled via 
Ising interaction. The field f(t) is varied with respect 
to time t and at time t = takes fo = /(0). The as- 
sociated controlled dynamical evolution is given by the 
Schrodinger equation 

j t \m) = -iH f (t)\m) (3) 

where we have assumed units such that H = 1. 

Our main objective is to prepare an Ising chain of 
length N with Hamiltonian (JlJ in one of the following 



III. SYMMETRIES, NON-CONTROLLABILITY, 
REACHABILITY 

The dynamics p| can be expressed in terms of the 
unitary process U(t) satisfying 

U(t) = -iH f {t)U{t) (6) 

with £7(0) = /. Denoting the solution of (|6| for a given 
control f(t) by U(t, f(t)), the reachable set of unitary op- 
erators 1Z is defined as the set of unitary matrices U(t, /) 
that can be generated by the dynamics Q in a finite time 
tf for some admissible control / £ J 7 , i.e., U £ 1Z if and 
only if there exists an admissible control f(t) such that 
U(tf,f(t)) — U. From control theory, we have [2D] : 
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Theorem 1. Let £ be the Lie algebra generated by 
span^ gjr { J[H + f(t)Hi]} = span{iJ , Hi}, called the 
dynamical Lie algebra. Then TZ = e . 

Hence, if £ = su(N) or £ = u(N) then we have 
K = SU(AT) or K = U(N) , the system is controllable 
in that we can implement any unitary operator U up to 
at most a global phase e 1 ^ . In this case any pure state 
is reachable from any other pure state, and more gener- 
ally, any two density operators with the same spectrum 
can be interconverted. Since the spectrum of any opera- 
tor is preserved under unitary evolution, this is the most 
we can hope for. In this sense, controllability is a suffi- 
cient condition for reachability. On the other hand, any 
system that possess symmetries will not be fully con- 
trollable. For the Hamitonian Hf(t), if there exists a 
Hermitian operator M such [M,H m ] = for m = 0, 1, 
then the Hamiltonians Hq and Hi are simultaneously 
block-diagonalizable. In this case the Hilbert space can 
be decomposed into orthogonal invariant subspaces T-Lk 
such that any initial state IV'(O)) € Tik remains in T-L^ 
under the evolution, regardless of what control f(t) we 
apply. Since no states outside Hi can be reached from an 
initial state in "Hi, decomposability immediately implies 
non-controllability. However, a target state \tl>d) may still 
be reachable from an initial state | -0(0)) if both belong to 
the same subspace. In particular, this is the case if the 
system is controllable on the relevant invariant subspace. 

Applying these results to our Ising chain subject 
to global control we see immediately that the sys- 
tem ([2]) possesses symmetries as both Hamiltonians Hq 
and Hi commute with the "A-parity" operator M — 
n i=1 X n . Hence, Hq and Hi are simultaneously block- 
diagonalizable. In our case M has two eigenspaces with 
eigenvalues ±1, spanned by 
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(e k ± e 2 N-fc+i) > , Vfc = 1, . 



. . (7) 

where is the basis vector with 1 in the fcth position 
and all other entries 0; changing from the basis {e^} to 
an eigenbasis of M simultaneously block-diagonalizes Hq 
and Hi . This shows that the Ising chain with global con- 
trol is not controllable and explains why additional re- 
sources are required to obtain universal QC [9HT3] . How- 
ever, there are many useful tasks that can be performed 
under the evolution ([3]). 

For N > 2 there are further symmetries and both 
blocks are further decomposable. There are various ap- 
proaches to decompose the Hilbert space into invariant 
subspaces that are not further decomposable. One ap- 
proach is to proceed as before and find the symmetry 
operators M on each subspace. If V is an eigenbasis 
of M then H m = H m V will be block-diagonal. This 
can be done recursively until no further symmetries are 
found for any of the blocks, leaving us with indecompos- 
able blocks. This approach becomes tedious, however, 
when there are many symmetries. Alternatively, we can 
calculate the eigenvectors of a linear combination of the 



Hamiltonians, H = aHo + /3Hi. Letting V be a unitary 
matrix whose columns are the normalized eigenvectors of 
H, we define an adjacency matrix A = (ay) with ay = 1 
if the absolute value of the (i, j)th element of the matrix 
V^HqV is greater than some threshold value S, and oth- 
erwise, and find the connected components of A, which 
define the respective invariant subspaces. The accuracy 
of this approach depends on suitable choice of a, /3 and S. 
Choosing S — 10~ 8 and a = 2, /3 = 3, we calculated the 
subspace decomposition for Ising chains up to N = 14. 

Having found a decomposition of the system into in- 
decomposable subspaces the next step is to verify if the 
initial and target states both belong to the same invariant 
subspace. For J < we verified numerically that both 
l^o") and the GHZ state \ipQ ') belong to the same invari- 
ant subspace for N — 2, ... ,14. To establish reachability 
the next step is usually to try to show that the system is 
controllable on this invariant subspace. For N = 2, this 
is easy. Changing the basis from {|0), |1)} to {|+), |— )} 
the Hamiltonians become 



Hn = 
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and similarly 

|Vtf) = [l,0,0,0] T , 



[1,0,0,1]VV2. 



belong to a two- 



This clearly shows that \ipo) and 
dimensional subspace Tl s spanned by the basis vectors e± 

(s) (s) 

and e4, on which we have H^ J = X and H{ ' = -2Z, 
showing that the dynamical Lie algebra generated is 
SU(2). Hence, the system is controllable on this sub- 
space and is reachable from \tPq). 

A similar approach allows us to establish reachabil- 
ity of the target state from the initial state for N = 3. 
One might therefore hope that the the system is con- 
trollable on the relevant invariant subspace T-L s for all 
N. Unfortunately, this is not true for higher dimen- 
sions. For N = 4, for instance, the smallest subspace 
H s that contains both \tpg) and \ipy ') and is invariant 
under the Hamiltonian has dimension 6, while the dy- 
namical Lie algebra £ generated by Hq and Hi on the 
entire space has only dimension dim£ = 16. This is 
strictly smaller than the dimension needed for controlla- 
bility on the subspace, which is 35 = dimsu(6) for full 
controllability and 21 = dimsp(3) for pure-state con- 
trollability for a six-dimensional subspace 1 . Hence, the 



1 Density operator controllability on a subspace of dimension N 
requires that we can generate the full Lie algebra of (trace-zero) 
Hermitian matrices u(JV) (su(A r )); if we are only interested in 
pure-state controllability it suffices if we can generate symplectic 
Lie algebra sp(2V/2). See [21] 
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system cannot be controllable on H s . Explicit calcula- 
tions for various N suggest that dynamical Lie algebra 
on the entire Hilbert space is a 2 -dimensional reducible 
representation of u(iV), and that the dynamical Lie alge- 
bra on the smallest invariant subspace that contains both 
|^|q") and \ip^} is an irreducible representation of su(N) 
or u(N). As for N > 3 the dimension of the subspace T-L s 
is greater than N , this implies non-controllability on T-L s 
for N > 3. Therefore, in general a different approach is 

needed to show that IV^'*) i s reachable from Simi- 
lar problems arise when trying to assess the reachability 
or non-reachability of the other GHZ states |Vv ) from 
H'o) or |Vo" >■ 

IV. CONTROL METHODS FOR GHZ 
GENERATION 

We now discuss three different methods for generat- 
ing GHZ states Q from the separable states \iJjq) or 
\iPq): (i) adiabatic passage; (ii) Lyapunov control and 
(iii) the optimal control theory. We shall see that adi- 
abatic passage demonstrates asymptotic reachability of 
certain GHZ states from certain product states. It is 
also has many benefits in that a simple field can achieve 
high fidelity and rather robust population transfer for 
spin chains of varying length. A general drawback of 
adiabatic schemes, however, is that the target state is ex- 
actly reachable only in the limit t — > +oo, and the time 
required to prepare the target state with sufficiently high 
fidelity can be long. This prompts the question whether 
we could do better using some form of optimal control 
design either Lyapunov control or global optimal control. 

A. Adiabatic Passage 

Assume we initialize the system in the ground state 

|^ /o) ) of H f (Q) = J[H + f Ht] for some /(0) = f . 
If we can show that a particular GHZ state in Q is the 
limit of the ground state of Ht(tj) as f(tf) —> 0, then it is 
possible to adiabatically transfer the system to one of the 
GHZ states Q. Indeed, such a result has been implic- 
itly shown for an Ising "chain" with periodic boundary 
conditions [22] using the Jordan- Wigner transformation, 
and we can easily show directly that this result is true 
for proper chains using the fact that Hq and Hi simul- 
taneously commute with the A-parity operator M, and 
thus that [Hf(t),M] = regardless of the choice of J 
and f(t). We clearly have 

M\^) = (+l) N \^+), (8a) 
M\^) = (-1) N \^). (8b) 

This shows that \iPq) always has positive parity, while 
|^>o") has positive parity for N even, and negative parity 
for N odd. The same must hold for the finite values of /o, 



i.e., the eigenstate |-0o ) nas positive parity if f < 0, or 
/o > and N even, and negative parity if fo > and N 
odd. As the parity is a conserved quantity, the adiabatic 

( f \ 

limit state \ip(tf)) must have the same parity as \ip ). 

It is easy to see that the intersection of the two-fold 
degenerate ground state manifold of JH with the +1 
( — 1) eigenspace of M is unique. For J < the ground 
state manifold of JHq is spanned by the GHZ states 
{{ip^) '■ k = 1,2}, and it is easy to see that IV^) 

(2) 

has positive, and \ip d '} negative parity. For J > the 
ground state manifold of JHq is spanned by the GHZ 
states {|^ fe ') : k = 3,4}, and it is easy to see that \ip^} 
has positive, and \ip d ) negative parity. Thus we have 

|^ X) ) £ Reach{|^ /o) ) : J < 0, J/ < 0, or 

J < 0, J f a > 0, N even}, (9a) 

|Vd 2) ) S Reach{|4 /o) ) : J < 0, Jf > 0, N odd}, (9b) 

|V^ 3) ) € Reach{|4 /o) ) : J > 0, J/ < 0, or 

J > 0, Jf > 0, N even}, (9c) 

|V^ 4) ) e Reach{|V^ /o) ) : J > 0, Jf > 0, N odd}. (9d) 

Although adiabatic passage provides a way to drive 
the system to a certain GHZ state \ipd) in Q, strictly 
speaking it only implies asymptotic reachability of \^pd) 
from the given initial state \ipQ ) as the Adiabatic Theo- 
rem provides that the error between the exact final state 
\ip(tf)} and \tpd) will go to zero only as tf — > +oo, i.e., 
that the fidelity 

F(t) = | (ip d \ip{t)) | 2 -> 1 as +oo, (10) 

if the field f(t) changes sufficiently slowly so that the 
rate of the change of the ground state energy ei , is small 
compared to the energy gap between the ground state 
and the first excited state, i.e., Ae = €2 — ei- This is 
not really a problem in practice as we are likely to be 
satisfied if we can get sufficiently close to the target state 
in a finite time, and similar adiabatic schemes have in- 
deed been proposed, e.g., for strings of neutral atoms 
in |24) . Moreover, under certain conditions asymptotic 
reachability implies finite-time reachability. 

Proposition 1. //e £ ; where £ is the dynamical Lie alge- 
bra, is compact then asymptotic reachability of\ipd) from 
an initial state IV'(O)) implies that \tpd) is reachable from 
IV'(O)) infinite time. 

Proof. If \ipd) is asymptotically reachable from 
|V>(0)) then there exists a path \ip(t)) such 
that lim t _> +00 \tp(t)) — \ipd)- We can choose a 
time sequence {t n } with t n — > +00 such that 
\ip(t n )) = U n \ij){Q)) \iji d ), where U n € e £ is a 
unitary process. Since {U n } is a sequence in the 
compact group e £ , there exists a converging subse- 
quence {U nk } such that U nk —> U G e £ , satisfying 
U\ip(0)) = IV'rf)- From Theorem [l] U S e £ , i.e., there 
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exists a dynamical trajectory during [0, T] for some finite 
T such that \ip{T)} = \ij) d ). □ 

According to [23J the dynamical Lie group is the di- 
rect product of an Abelian Lie group and a semi-simple 
compact Lie group. Numerical computations of the gen- 
erators of the Abelian Lie group for various N suggest 
that in our case, the Abelian group is compact, and thus 
that e £ generated by Hf{t) is compact, and hence we 
have finite-time reachability. 

To assess the performance of adiabatic passage we turn 
to simulations. Although the choice of f(t) does not mat- 
ter in theory, provided it varies sufficiently slowly and 
vanishes as t — > +oo, in practice we are usually interested 
in preparing a sufficiently close approximation to the tar- 
get state in as little time as possible, and in this case the 
choice of f(t) does matter as a comparison of two simple 
controls, a linearly decreasing field f(t) = fo(l — t/tj) 
for t G [0, t/], and a decaying exponential, f(t) — /oe~ pt , 
in Fig. [T] shows. The results of these simulations suggest 
that the latter choice is preferable in terms of speed and 
robustness. This can partly be explained by comparing 
E(t) = ei/Ae. For the linear field E(t) is negligible for 
most of the pulse duration and spikes towards the end 
of the pulse, mirroring the sharp increase in the popu- 
lation of the target state near the target time (Fig. [I] 
left). For the exponentially decaying field f(t) drops 
to / ss 1 much faster and spends more time in the re- 
gion < / < 1, where most of the interesting evolu- 
tion takes place. Also, for the exponential field we have 
f(t)/f(t) — p, i.e., E(t) s» —2.5p is approximately con- 
stant until the field has dropped to / « 1, while for the 
linear field f(t)/f(i) = p/f(t), where p is the slope, and 
thus E{t) « -2.5/Lt//(t) will be negligibly small for f(t) 
large. 

In both cases the choice of /o is dictated by practical 
concerns. Fig. [2] shows the ground state energy e\ and 
the energy gap Ae for Ht(t) with J < 0, as a function 
of field strength /. It shows that for / > 1 the energy 
gap increases with /. Assuming the rate of cooling to 
be proportional to the energy gap, this compels us to 
choose f as large as possible to ensure efficient cooling 
and initial state preparation. It also ensures that the ini- 
tial state is close to the desired initial state l^o") as the 

error 1 — |(V'o"l'0o^ O ^)| 2 decreases quadratically in / . At 
the same time, to maintain adiabaticity, the field must 
decrease slowly, and thus the time required for the field 
to decay to a certain value close to zero increases propor- 
tionally. For fixed /o the asymptotic value of the fidelity 
depends on the decay rate p. In general it decreases as p 
increases. At the same time, the time required to reach 
a certain target fidelity (below the asymptotic value) de- 
creases with increasing p. Thus, if we are only interested 
in achieving a certain target fidelity of say 99%, there will 
be an optimal value of the decay rate p that achieves 99% 
transfer the shortest amount of time. Fig. [3] shows the 
log-error log 10 (l — F) for chains for different length N for 
an exponentially decaying field with p = 0.1. We note 



that the time when the fidelity crosses the 99% thresh- 
old remains in a narrow range of [44, 49] for a significant 
range of N. The onset of oscillations in the evolution 
of the fidelity (population of the target state) for larger 
N suggests that non-adiabatic effects arise for this p, re- 
sulting in population transfer from the (instantaneous) 
ground state \ipi(t)) of H(t) to excited states. This is 
confirmed by the population plot 1 — pi(t) for N = 11 
(Fig. [3| inset), and shows that we must reduce the decay 
rate p for longer chains to maintain adiabatic evolution. 



B. Lyapunov control 

A simple way to solve certain optimal control prob- 
lems is construct a Lyapunov function for the dynamics 
([3]). A natural candidate for a Lyapunov function is a 
monotonic function of the Hilbert-Schmidt distance for 
density operators such as 



v{p ) p d ) = \\\p-p d r = \Tr[{p- Pd n 



(ii) 



satisfying V > which equality holds if and only if 
Pd = p- The essential idea of Lyapunov control is to de- 
sign appropriate control dynamics such that V becomes 
a Lyapunov function, i.e., V keeps decreasing along every 
trajectory of p(t). This can be realized by choosing 

/(*) - f(p(t),pd) = K'&a^i.pdH*))- (i 2 ) 

Then for V(t) = V(p(t),p d ), we have 

V(t) = -f(t) Tr([iHi,p d ]p(fi) = - K f(t) 2 < °> (I 3 ) 

i.e., p(t) evolves towards p d . Ideally, if V(t) — > as 
t — > +00, we have p(t) — > p d , but this does not always 
hold in general. However, the LaSalle invariance princi- 
ple [25] ensures that every solution p(t) under ( 12 ) con- 
verges to a set, often known as the LaSalle invariant set, 
and it has been shown that the control design above ren- 
ders the target state p d almost globally attractive if (i) 
Ho is strongly regular and (ii) Hi is fully connected [26j . 
Condition (i) requires Hq has distinct transition frequen- 
cies between any pair of energy levels in the smallest 
invariant subspace H s containing po and p d . 

For N = 2 spins, with the initial and the target states 
chosen as before, it can be verified that the two con- 
ditions above hold for the dynamics restricted to T-L Sl 
and the GHZ state p d is almost globally attractive on 
W. s . Hence the Lyapunov control pulse (12) will steer 
the system from \tp^ ) to p d . Moreover, as demonstrated 
in |26j . the Lyapunov control design ensures that p(t) 
converges "exponentially" to p d , and is insensitive to tim- 
ing errors — unlike geometric control schemes that require 
precise on-off switching of the control fields and ampli- 
tude control, for instance. Unfortunately, for longer Ising 
chains, convergence of p(t) to p d is no longer assured, as 
the energy levels of H on H s are equally spaced, i.e., 
condition (i) does not hold, and for N > 4 full connec- 
tivity is lost. Theoretical analysis [57] shows that in this 
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FIG. 1: (Color online) Adiabatic control (top) and corresponding relative energy gap E(t) = ei(i)/Ae(t) and fidelity F(t), for 
a linearly varying field, f(t) = 10(1 — t/100), left, and and exponentially varying control, fit) = 10e~ 05 ', (right). The fidelity 

Fit) corresponding to the projection onto the GHZ state in both cases asymptotically approaches a limiting value of 

« 1, but for the linear control it increases sharply near the end of the pulse. For the exponentially decaying control field the 
increase begins much sooner, is more gradual, and higher fidelities can be achieved in shorter time. 
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FIG. 2: (Color online) Ground state energy e 1 (solid line) 
and energy gap (dashed) Ae = £2 — ei for a chain of length 
10 with Hamiltonian Hf, given in Eq. l|l]), and J < as a 
function of the applied field /, for the dynamics restricted to 
the invariant subspace H s containing the GHZ state IV'd }• 
For large fields the Hi term dominates and the ground state 
energy and energy gap depend approximately linearly on the 
field: ei(/) « -10/, Ae(/) ^4/. For < / < 2 the inset 
suggests an approximately quadratic dependence of both the 
ground state energy and energy gap on the applied field. 



case the invariant set is large, and most solutions con- 
verge to "limit cycles" a finite distance from the target 
state pd- Nonetheless, simulations suggest that we can 



ensure p(t) with p(0) = IV'o /(V'o I converges to a point 
very close to pd by carefully tuning the so-called feed- 
back strength k in (12), as illustrated in Fig. [4] We find 
that (i) for a given N the final fidelity achieved usually 
increases with k; (ii) for a given k the final fidelity de- 
creases with increasing N. This can be explained by the 



FIG. 3: (Color online) The log-error, i.e., log of population 
not transferred to the target state, as function of time for 
Ising chains (Jl]) with J < of length N = 2, , . . 11 starting in 
|V>o /o) ), shows that the field f(t) = 10e" o lt achieves 99% pop- 
ulation transfer to the target state for chains of length N = 2 
to 10 in a narrow time window to. 99 G [44,49]. For larger 
N oscillations indicate non-adiabatic effects, which prevent 
reaching 99% fidelity for N = 11 with this field. The popu- 
lation of the upper levels 1 — pi(t) (inset) clearly shows that 
the pulse induces some population transfer to these states, 
occuring, as expected, following the peak in E(t), indicating 
that the peak value is too large to fully maintain adiabaticity. 



fact that the dimension of the center manifold surround- 
ing pd increases rapidly with ^V. Moreover, to achieve 
higher fidelities for larger k, the control time increases 
rapidly as well. Therefore, considering the transfer time 
required to achieve a certain target fidelity close to 1, 
Lyapunov control has a significant edge for small N, but 
the advantage disappears for longer chains (see Fig. |4| . 
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FIG. 4: (Color online) Examples of Lyapunov controls (solid 
lines) and "distance" V(o(t), pd) (dashed) from the GHZ state 
IV'd ) f° r Ising chains |l| with J < 0, of different length N, 
assuming the system is initially prepared in the ground state 
\ipQ ) with /o = 10. The Lyapunov control pulses achieve the 
desired state transfer in much shorter time than the adiabatic 
controls for small N, and the approximately linear decrease 
of the distance on the logarithmic scale, suggests that the 
distance decreases exponentially. For longer chains, however, 
the feedback strength k has to be increased and the effects of 
"attractive" limit cycles interfere with convergence, resulting 
in control pulses comparable or even longer than the corre- 
sponding adiabatic transfer pulses. 

Thus, the Lyapunov control design is only effective for 
short Ising chains. 

C. Optimal Control 

Lyapunov control can be considered as a kind of op- 
timal control in that the distance from the target state 
is monotonically decreasing with time. This form of op- 
timal control is sometimes referred to as local optimal 
control because at each point in time the control design 
is based only on information about the current state of 
the system and the target state. From a computational 
point its main advantage is that the control pulses can be 
calculated directly in a non-iterative fashion, but as the 
last section shows, for more complex problems such as 
longer chains this approach is not sufficiently powerful. 
An alternative method is to take a global approach, spec- 
ify a target time tf, and attempt to maximize the fidelity 



F(tf) by globally varying the control pulse f(t) in the al- 
lowed control function space. In practice, optimal control 
problems over function spaces can generally be solved 
only numerically, by parameterizing or discretizing the 
control f(t). The simplest and most common approach 
is to subdivide the time interval [0,i/] and approximate 
the control f(t) by a constant on each subinterval 
This results in an optimization problem over R^, where 
K is the number of time intervals, which can be solved, 
e.g., by starting with an initial trial field f^°\t) and iter- 
atively refining f^ n \t) such that the fidelity at the target 
time F^ n '{tf) monotonically increases as a function of the 
iteration n. The crucial part of this procedure is the way 
f( n '(t) is updated in each iteration. We adopted here 
a quasi-Newton method developed by Boyden, Fletcher, 
Goldfarb, and Shanno [55] • A more detailed discussion 
of the optimization process can be found in [29] , 

The adiabatic and the Lyapunov control pulses found 
earlier provide upper bounds on tf to reach a certain 
final fidelity F < 1, and our previous reachability consid- 
erations for Ising chains with Hamiltonian Ht(t) suggest 
that there exists a solution fit") with F(tf)) = 1 for a 
finite time t f , although the proof does not give any hint 
as to the best driving field f(t) or how to discretize the 
control pulse. Although it can be shown, roughly, that if 
the target state is reachable by some admissible control 
field then there also exists a piecewise constant control 
that achieves the same task, the theorem again does not 
tell us how many time steps are needed. If the resolution 
is too low, i.e., At = tf/K is too large, the target fidelity 
may not be achievable. An interesting question for opti- 
mal control therefore is how fast we can hope to achieve 
transfer to the target state, and what time resolution of 
the field is required. Can we do better than adiabatic 
control? 

Before we address these questions a final issue that 
needs to be considered is constraints. For example, the 
amplitude of the control pulses will usually be limited by 
what can be achieved experimentally. Whether to impose 
constraints and what kind depend on the details about 
the system and the implementation. For typical values in 
NMR experiments J is usually a few hundred Hz, while 
we can easily achieve fields up to 50kHz for liquid-state 
NMR and is a few hundred kHz for solid-state NMR [5], 
which would suggest a reasonable upper bound on the 
field amplitudes of perhaps 10 3 | J|. Simulations for this 
problem suggest that bounds of this magnitude can usu- 
ally be neglected as the unconstrained optimization solu- 
tions almost always satisfy these constraints, and uncon- 
strained optimization is computationally more efficient 
than constrained optimization. 

Using an unconstrained optimization based on the 
above-mentioned quasi-Newton method, we calculated 
the optimal pulses for Ising chains of length N = 2 to 
12. Assuming J < and starting with the initial state 
\4>q), we calculated optimal controls for different values 
of the target time t f and time steps K, choosing the GHZ 
state \i>} ) as a target state. For chains up to N = 12 




Length of chain N 



FIG. 5: (Color online) Lower bounds on the time tf (solid line; 
y-axis linear scale) required to achieve 99% transfer fidelity as 
a function of chain length for TV = 2 to TV = 12, derived from 
optimal control simulations. Note the linear dependence oft/ 
on TV, despite an exponential increase in the dimension of the 
relevant subspace H s (dashed line; y-axis logarithmic scale). 



we were able to find controls that achieve at least 99% 
fidelity for tf(N) — 0.65TV and K = 3TV, and in general 

only one or two runs with either /> ' = 1 or a random 
number in [0, 1] for k = 1, . . . , K were needed for the op- 
timization to succeed in finding a suitable control. This 
shows that the transfer times for the optimal controls 
are much shorter than those for the adiabatic pulses - 
by about one order of magnitude - and the transfer time 
appears to be increase linearly with the chain length TV, 
at least up to TV = 12 as shown in Fig. [5j This is quite 
surprising when one considers that the dimension of the 
smallest invariant subspace H s that contains the GHZ 
state \ipj ) increases exponentially in TV. A possible ex- 
planation for this linear dependence lies in the dimension 
of the reachable set. The dynamical Lie algebra on the 
subspace T-L s appears to be a high-dimensional represen- 
tation of su(TV) or u(TV), suggesting that the reachable 
set starting with a pure initial state is the homogeneous 
space U(TV)/[U(1) (g> U(TV — 1)], which has dimension 
TV 2 - 1 - (TV - l) 2 = 2TV - 2 30,. This suggests that the 
reachable set is a (2TV — 2)-dimensional manifold, which 
may explain the apparently linear dependence in the ob- 
served bounds on the transfer times despite the expo- 
nential increase in the dimension of the subspace ~H S it is 
embedded in (Fig. [5]). 

Fig. [6] shows an example of a typical optimal control 
field fjt) for a chain of length TV — 10 and the corre- 
sponding evolution of the system. The pulse appears 
well-behaved and quite feasible both in time and fre- 
quency domain. The corresponding evolution of the sys- 
tem shows that the fidelity does not increase monoton- 
ically, as is the case for an ideal adiabatic or Lyapunov 
control pulse. The figure also shows the evolution of the 




1 2 3 4 5 



time (in units of 1/|J|) 

FIG. 6: (Color online) Optimal control pulse (top) and evolu- 
tion of the ground state population Si (or fidelity F(t)) and 
other excited states population (Sa, S3, Sio) for an Ising 
chain of length TV = 10 with J < 0, starting in the state 
|-0q ). The fidelity corresponds to the projection onto the tar- 
get GHZ state \ip^ ). Inset shows spectrum of optimal pulse. 



population of the eigenspace of H . For TV = 10, H , re- 
stricted to the relevant subspace W s of dimension 252, has 
10 distinct eigenvalues with corresponding eigenspaces 
S n of varying dimensions. The population of the ID 
ground state manifold Si corresponds to the fidelity F(t). 
Analysis of the population evolution shows that the in- 
crease in the fidelity F(t), or the ground state popula- 
tion Si, is preceded by an increase in the populations of 
the first and second excited state manifold, and these in- 
creases are in turn preceded by peaks of the populations 
of the eigenspaces S4, S5 and S&, in this order. This be- 
havior can be explained in terms of the subspace coupling 
induced by the interaction Hamiltonian Hi . For the dy- 
namics restricted to the subspace T-L s , Hi couples each 
eigenspace of Hq only with its first and second neighbor, 
i.e., the ground state manifold is directly coupled only to 
the first and second excited state manifold, and so forth. 
Therefore population cannot be directly transferred from 
e.g., 1S10 to Si, but must pass through several intermedi- 
ate levels. Note that the maximum number of intermedi- 
ate levels to be traversed is linear in TV as the restriction 
of Hq to the subspace T-L s containing the GHZ state \ipd ) 
for a chain of length TV has TV distinct eigenvalues. If the 
system starts in the state \ipo) or a state close to it, all 
of these eigenspaces are initially populated. This does 
not fully explain the behavior of the optimized dynam- 
ics as there are many different excitation pathways and 
optimal control attempts to maximize constructive inter- 
ference between all different paths leading to the desired 
outcome, but it suggests an alternative explanation for 
the apparently linear dependence of the minimum trans- 
fer time on TV, at least for chains up to length TV = 12, 
despite the exponential increase in the dimension of T-L s 
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FIG. 7: (Color online) Final fidelity achieved for control pulse 
f(t), which is optimal for J < and the initial state | V^o" ^ ; 

while the actual initial state is the ground state | "0o^° ^ ) ' as a 
function of fo. The error between the ideal (dashed line) and 
real fidelity (solid line) goes to zero as fo increases. 



FIG. 8: (Color online) The final f idel ity achieved by optimal 
control for a thermal initial state ( 14 1 instead of a pure state 



decreases with temperature, with a faster drop-off for longer 
chains, as theoretically predicted, but high transfer fidelities 
can still be achieved for reasonably low temperatures. 



from 2 for N = 2 to 924 for N = 12. 



V. EFFECT OF IMPERFECTIONS 

So far it was assumed that we have an ideal chain with 
uniform couplings and can perfectly initialize the system 
in a pure separable state such as \ipQ ) or the ground state 

of I^q^) for a fixed large value of fo by cooling, and sys- 
tem inhomogeneity and the effect of interactions with an 
environment were neglected. To assess the robustness 
of various control strategies with regard to such imper- 
fections in the context of this problem we shall consider 
how the effectiveness of the control fields they produce 
for ideal systems is diminished by various imperfections. 



A. Imperfect Initialization 

Preparing the system in the ideal state | V^o" ) is n °t 
always an easy task since we need to achieve J fo — > — oo. 
We have already discussed using a finite fo for adiabatic 
passage and it was shown that provided fo is sufficiently 
large, we always can realize one of the GHZ states Q. 
For Lyapunov and optimal control strategies we may have 
been told that the system is initialized in the state | V^o" ^ ' 
while it was really initialized in the ground state IV'o ) 
for some finite value of fo- In this case f(t) is optimal 
for the initial state l^o") and using this field for the real 

initial state \ipQ ) does not give us the maximal fidelity. 
The effect of such errors is easy to analyze. Assume we 
have an ideal (optimal, adiabatic) control f{t) that gives 
rise to an evolution U{t) such that C(t/)|^) = \4>a ) 



for some finite tf or tf — > +oo. Let IIo = |V'o~)('0o"| be 
the projector onto IV'o")- Then we have 

= u \4 M ) + nM fo) ) = + ^M fo) ) 

with Co = (V'dlV'o^ '')- Furthermore 

+ (^\U(tM\4 M )=co 

as we have U(tf)\ipg) = IV^); an d U(tf) is unitary. So, 
the maximum transfer fidelity is simply given by |co| 2 , 
the overlap of the actual initial state with the assumed 
initial state, as illustrated in Fig. [7j For adiabatic con- 
trol any deviation of the initial state from the ground of 
the Hamiltonian at time t — will limit the maximum 
transfer fidelity with the bound given by 1 — | Co | 2 - For 
optimal control the error can be decreased by determin- 
ing the initial state more accurately, e.g., using system 
identification techniques such as state tomography, and 
incorporating this information in the optimization. 

B. Thermal Effects 

Assuming the system is initialized by cooling in the 
presence of a strong field in the x-direction, another 
source of error are thermal effects, in particular the fact 
that zero temperature is not practically achievable. In 
general finite temperature effects result in an initial state 
that is a thermal mixed state 

e -H f (0)/k B T 

P(°) - Tl[e - Hm/kB Ty ( 14 ) 
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where T is temperature and is the Boltzmann con- 
stant. As the populations of the eigenstates cannot 
change under adiabatic passage, wc sec immediately that 
the maximum projection onto the target GHZ state we 
can achieve is given by the ground state population 

Wl = <V>o C/o V(0)l4 /o) >- (15) 

The initial population of the ground state depends both 
on the temperature T and the energy gap between the 
ground state and the excited states of i?/(0). As the 
energy gap Ae increases roughly linearly with the field 
strength / , as we have seen earlier, this explains why it 
is desirable to cool in the presence of a strong field. 

Based on the results in [29^ we might assume that the 
transfer fidelity for optimal control could be improved by 
starting with thermal initial state rather than the ground 
state. Unfortunately, this is not the case here because 
unlike in [29] the observable we are optimizing is a rank- 

1 projector onto a pure state, \ipd){ipd \, i.e., it has a 
single eigenvalue of 1 and all other eigenvalues are zero. 
Therefore, the transfer fidelity is bounded above by [3T] 

F{tf) = {^ ) \p(t f )\^ ) )<w 1 . (16) 

Fig. [H] shows the actual fidelities achieved by optimal 
control for Ising chains of length N = 6, 7, 8 with J < 
as a function of temperature T, assuming the system is 
initialized in a thermal ensemble of Hf(0) with fo = 10. 
The results are in line with the expected decrease of the 
ground state population w\ as a function of the initial 
temperature T, and the fact that w\ decreases faster for 
longer chains. In Fig. [9j for Ising chain with N = 6 
and J < 0, assuming initialized in the thermal mixed 
state, we have plotted how the final fidelity changes with 
respect to temperature for the three different methods. 
The relative flatness of the curve for low temperatures 
suggests that all of these control methods can achieve 
high-fidelity state transfer for reasonably long chains and 
sufficiently low temperatures however, but the optimal 
control pulses appear more robust with regard to thermal 
fluctuations than the two other methods. 



C. Disordered Chains 

In reality it is impossible to make an absolutely uni- 
form chain. There are always fluctuations in the spin 
coupling which make the chain disordered. To assess 
the effect of such inhomogeneity we consider the case 
where the couplings between neighboring sites are ran- 
domly perturbed around the average value J 

N-l N 

H' f {t) = ^J(1 + S n )Z n Z n+1 + f(t) x n, (17) 

n— 1 n— 1 

where 5 n G [— a. +a] is a random variable with uniform 
distribution around zero. Since the random couplings are 
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FIG. 9: (Color online) Final fidelities versus temperature for 
different control methods for a chain of length N = 6: (i) op- 
timal control, (ii) adiabatic control, and (iii) Lyapunov con- 
trol, when the control pulses generated assuming the system 
is initially in the ground state of Ho + foHx, are applied to 
initial states that are thermal ensembles of different temper- 
ature. The decrease in the final fidelity as a function of the 
temperature is the similar for all three methods. 



not known we can only study their average effect over en- 
tanglement generation. Thus, we consider the evolution 
of the initial state \4>q) under action on the perturbed 
H'f{t) for a pulse f it) optimized assuming the Hamilto- 
nian Hf(t). We repeat the experiment over 100 random 
perturbations and then take the average value of the final 
fidelity over all results. Fig. [10] shows the average fidelity 
as a function of the parameter a of the fluctuations in the 
coupling strength for different methodologies. The figure 
shows that adiabatic passage is very robust against dis- 
order in the chain. This is to be expected as the ground 
state of the Hamiltonian H'Jt) is very similar to Hf(t) 
for all times when the fluctuations are small, and there- 
fore adiabatic passage is always able to steer the Ising 
chain from to the final GHZ state IV^) indepen- 
dent of the exact choice of the pulse fit). Lyapunov and 
optimal control pulses on the other hand rely on dynamic 
and interference effects to achieve faster transfer, and the 
temporal shape of the optimal pulse is thus more strongly 
dependent on the form of the Hamiltonian including the 
coupling strength. Thus, the optimal control pulses are 
more susceptible to disorder although the optimal control 
pulses appear to outperform adiabatic control for fluctu- 
ations up to a w 0.03. It is also interesting to note that 
the more efficient control pulses obtained from global op- 
timization techniques are more susceptible than their less 
effective Lyapunov control cousins. Again, the perfor- 
mance of optimal control schemes can be significantly 
improved if the actual couplings can be estimated more 
accurately using system identification techniques [32 or 
closed-loop adaptive strategies [33]. Alternatively, when 
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FIG. 10: (Color online) Comparison between different control 
methods for GHZ generation for a chain of N = 6 in terms 
of a, the strength of disorder. The target time was tf = 
50J" 1 for adiabatic passage and tf = 15J -1 for Lyapunov 
and optimal control. 



FIG. 11: (Color online) Effect of dephasing on different meth- 
ods for GHZ generation for N = 6. The longer pulses required 
for adiabatic passage leave it more susceptible to dephasing. 
tf — 50 J^ 1 for adiabatic passage, tf = 15 J -1 for Lyapunov 
and optimal control. 



accurate estimation of the couplings is infeasible or im- 
possible, e.g., because we wish to control an ensemble of 
systems with slightly different couplings, the performance 
of optimal pulses can (sometimes) be improved by opti- 
mizing ensemble averages for a collection of systems as in 
[34] . although this may not be feasible for systems with 
more than a few spins as optimizing over a collection of 
systems is computationally expensive, with each function 
evaluation requiring the numerical solution of the time- 
dependent Schrodinger equation for many systems. 



D. Decoherence Effect 

Another major problem in practice is that it is impos- 
sible to isolate the system from its environment. Any 
interaction with environment tends to disturb the evolu- 
tion of the system. The precise effect of the environment 
clearly depends on the type of interaction. We shall as- 
sume that the system-environment interaction is weak 
and Markovian and modelled by a Lindblad equation 

p(t) = -i[H f (t),p(t)]+C(p(t)), (18) 

where C(p) corresponds to dissipative effects and p(t) is 
the density matrix of the system. We shall focus here on 
dephasing, which destroys the coherence of the system, 
as it is common for many physical systems and the typi- 
cal dephasing times for most physical systems are much 
shorter than other relaxation rates. 

To assess the effect of decoherence we solve the Lind- 
blad equation using the appropriate adiabatic or opti- 
mized pulse fit) obtained for the noiseless system. For a 



dephasing noise we consider the specific model 

JV 

£(P(*)) = -7X>W - Z np(t)Z n }, (19) 
n=l 

which corresponds to localized dephasing of individual 
spins. Fig. [TT] shows the fidelity achieved as a function of 
noise strength 7 for all methods. As expected, the fidelity 
decays as the noise strength increases for all strategies 
but adiabatic passage is substantially more susceptible 
to the dephasing noise of the form considered. This is 
mainly due to the times involved: the transfer times for 
the adiabatic pulses tend to be about one order of mag- 
nitude greater than those for the optimal pulses, giving 
dephasing more time to act and destroy the coherence. 
Optimal control is more robust, and theoretically, it may 
be possible to improve the performance of optimal con- 
trol by taking decoherence effects into account in the op- 
timization progress, although this is computationally de- 
manding as it require full density matrix optimzation for 
a density matrix of Hilbert space of dimension 2 , as 
opposed to pure-state optimization on a subspace, which 
can be done far more efficiently. 



VI. CONCLUDING DISCUSSION 

We have studied Ising chains subject to a single a 
global control such as a magnetic field in the x-direction 
with regard to preparing them in an entangled GHZ 
state. Due to the existence of multiple symmetries the 
system is not controllable, and the Hilbert space decom- 
poses into subspaces of varying dimensions that are in- 
variant under the dynamics. Our analysis further shows 
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that even on the invariant subspaces the dynamics is not 
controllable for all but the smallest subspaces for TV > 3, 
and that the GHZ state of interest always lies in one of 
the largest invariant subspaces on which the dynamics 
is not controllable except for N = 2,3. Nonetheless the 
GHZ state is (asymptotically) reachable from an easy- 
to-prepare initial state by adiabatic passage, and under 
certain conditions finite-time reachability from a product 
state such as | + . . . +) can be inferred for chains of arbi- 
trary length TV. Motivated by this positive reachability 
result, three different control strategies — adiabatic con- 
trol, Lyapunov control and optimal control — were consid- 
ered to prepare a GHZ state, and their advantages and 
disadvantages discussed, serving as a reference for future 
experimental implementations. 

Each method has its own advantages and disadvan- 
tages. Broadly speaking Lyapunov control is a simple 
form of optimal control design, which produces pulses 
that are simple and quite effective for short chains but the 
method struggles for longer chains, and for TV > 4 both 
adiabatic and global optimal control seem to be generally 
superior. In terms of the total time required to prepare 
a GHZ state with high fidelity for a given chain of length 
TV, the control pulses found by global optimal control 
proved the most efficient, with the time required being 
up to an order of magnitude less than the time-scale for 
adiabatic transfer. The short pulse durations also con- 
fer greater robustness in the presence of local spin de- 
phasing compared to the adiabatic transfer scheme, and 
the optimal control pulses appear slightly more robust 
with regard to finite temperature effects. The adiabatic 
transfer scheme, on the other hand, is more robust with 
regard to system inhomogeneity such as unknown fluc- 
tuations in the J-coupling between spins. This is to be 
expected considering that small system perturbations do 
not change the ground state of the system appreciably, 
while such perturbations can alter the interference be- 
tween different excitation pathways that optimal control 
designs aim to exploit. The performance of the latter 
can be improved by system identification or closed-loop 
adaptive strategies, and in some cases by explicitly taking 
inhomogeneity or decoherence into account in the opti- 



mization. 

Both the adiabatic pulse and the Lyapunov pulse have 
analytic expressions, and especially for the adiabatic con- 
trol, we can choose the form of the slowly varying pulse 
without solving the dynamical equation. It is interesting 
in this regard to note that although a linearly decreas- 
ing field would appear to be the simplest choice for adi- 
abatic control, simulations and analysis suggest that a 
field of the form /oe _M * is generally preferable, resulting 
in both faster transfer and increased robustness. While 
the pulse shapes of the optimal pulses are rather random 
and more complex than the corresponding adiabatic con- 
trol fields, the amplitude range and spectral bandwidth 
of the pulses are quite narrow, and certainly appear to 
be within experimentally accessible limits. One major 
drawback of optimal control is that, unlike for the adi- 
abatic and Lyapunov control designs, there are no ex- 
plicit expressions for the optimal pulse. Instead the pulse 
must be computed by numerically, and the complexity 
and computational overhead of calculating the optimal 
pulses increases rapidly with the length of the chain. Al- 
though the dynamics can be restricted to a subspacc, un- 
fortunately, the GHZ state of interest belongs to one of 
the largest invariant subspaces H s , whose dimension in- 
creases exponentially in TV with dim'H s w 2^/4 for large 
TV, which substantially increases the computational com- 
plexity. Surprisingly, despite this exponential increase of 
the subspace dimension, the minimum time required to 
prepare a GHZ state given an initial product state such 
as | + . . . +) using optimal control appears to be actually 
linear in TV. Such fast production of a multi-spin GHZ 
(or Cat or NOON) state using a minimal global field on 
an Ising chain will be highly valuable in the fields of en- 
hanced sensing and quantum information technology. 
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